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Abstract 

The LHCb Vertex Locator (VELO) is a shicon strip detector designed to reconstruct 
charged particle trajectories and vertices produced at the LHCb interaction region. 
During the first two years of data coUection, the 84 VELO sensors have been exposed 
to a range of fluences up to a maximum value of approximately 45 x 10^^ 1 MeV 
neutron equivalent (IMeVueq). At the operational sensor temperature of approx- 
imately — 7°C, the average rate of sensor current increase is 18 |xA per fb~^, in 
excellent agreement with predictions. The silicon effective bandgap has been deter- 
mined using current versus temperature scan data after irradiation, with an average 
value of Eg = 1.16ib0.03ib0.04eV obtained. The first observation of n'''-on-n sensor 
type inversion at the LHC has been made, occurring at a fluence of around 15 x 10^^ 
of IMeVucq. The only n~^-on-p sensors in use at the LHC have also been stud- 
ied. With an initial fluence of approximately 3 x lO^'^ 1 MeVugq, a decrease in the 
Effective Depletion Voltage (EDV) of around 25 V is observed, attributed to oxy- 
gen induced removal of boron interstitial sites. Following this initial decrease, the 
EDV increases at a comparable rate to the type inverted n''~-on-n type sensors, with 
rates of (1.43 ±0.16) x 10"^^ y/ i MeVueq and (1.35 ±0.25) x 10"^^ y/ ^ ^gVueq 
measured for n~^-on-p and n"''-on-n type sensors, respectively. A reduction in the 
charge collection efficiency due to an unexpected effect involving the second metal 
layer readout lines is observed. 
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Figure 1: A schematic representation of an R-type and a (p-type sensor, with the routing hues 
orientated perpendicular and parallel to the silicon strips, respectively. 

1 Introduction 

The VErtex LOcator (VELO) is a silicon strip detector positioned around the proton- 
proton interaction region at the LHCb [1] experiment. To obtain the precision vertexing 
required for heavy-flavour physics, the closest active silicon sensor region is located 8.2 mm 
from the beam axis, while the silicon edge is located at a distance of 7 mm. For the lumi- 
nosity delivered by the LHC in 2010 and 2011, the VELO was exposed to higher particle 
fiuences than any other silicon detector at the LHC. Careful monitoring of radiation dam- 
age to the sensors is essential to ensure the quality of data for LHCb physics analyses and 
to provide information relevant to the eventual detector replacement and upgrade. 

During proton injection and energy ramping the LHC beams are wider and less stable 
than the beams used for data taking. To prevent damage to the silicon sensors, the VELO 
consists of two halves retractable by 29 mm in the horizontal plane. Each half contains 42 
half-disc shaped silicon-strip sensors. When the beams are in a stable orbit, the two VELO 
halves are closed such that the colliding beams are surrounded by the silicon sensors. Half 
of the sensors have strips orientated in an approximately radial direction (0-type) and 
the other half perpendicular to this (R-type), as shown in figure [l] A detector module 
consists of an R-type and a (p-tjpe sensor glued to a common support in a back-to-back 
configuration. Track coordinates are measured using the charge collected by the two 
sensors in a module. All but two of the VELO sensors are oxygenated n"'"-on-n sensors, 
consisting of an n-type implant on a n-type bulk with a backplane p"'"-type implant. Two 
oxygenated n~^-on-p silicon sensors are installed at one end of the VELO, intended to be a 
test of one of the leading LHC silicon- upgrade candidates in an operational environment. 
A summary of the silicon sensor properties is given in table [Ij 

Each implant is read out via a capacitively coupled first metal layer running along 
its length. The R-type strips and inner 0-type strips do not extend to the outer region 
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of the sensor. Therefore each strip is connected via a metal routing hne to the edge of 
the sensor where the readout electronics are located. The routing lines are referred to as 
the second metal layer and are insulated from the bulk silicon and first metal layer by 
3.8 ± 0.3 |a.m of Si02. For R-type sensors, the routing lines are positioned perpendicular 
to the sensor strips, whilst for the 0-type they are positioned directly above and parallel 
to the strips. 

This paper presents studies of radiation damage in the VELO sensors using data col- 
lected during 2010 and 2011. Section[2] describes studies of sensor currents as a function of 
voltage and temperature. In section [3] the effects of radiation are monitored by measuring 
the charge collection efficiency and noise as a function of bias voltage. An unexpected de- 
crease in clustering efficiency due to an effect involving the second metal layer is described 
in section HI The results from the various studies are summarised in section |5l 



2 Current evolution 

The leakage current in a silicon sensor varies linearly with fluence for a wide range of silicon 
substrates [2] . This predictability provides a simple and accurate method of relating sensor 
currents to the amount of particle fluence accumulated by a sensor. This section presents 
an analysis of sensor currents in order to monitor radiation damage in the VELO. 

2.1 Interpretation of sensor current measurements 

Each VELO sensor is biased through a high-voltage connection to the backplane, and a 
common ground connection to the bias line and innermost guard ring on the strip side. 
Therefore a single current is measured, corresponding to the current drawn by the entire 
sensor. The raw measured currents are shown as a function of time in figure [2| for sensors 
operated at the nominal bias voltage of 150 V and at a mean temperature of approximately 
— 7°C. During the early running period, the majority of the sensors had small currents, 
with a few exceptions that had currents of up to 40 f^A. With sensor irradiation the 



Table 1: The VELO sensor design parameters. The sensor position along the beam axis is given 
relative to the beam interaction region. 



Parameter 



Value 



Silicon thickness 

Strip pitch 

Strip width 

Routing line width 

Inner silicon edge to beam axis 

Radial distance of active strips from beam axis 

Sensor position along beam-axis 

Oxygen enhancement 



300 i^m 
40-120 M-m 
11-38 M-m 
~11 |j,m 
7 mm 
8.2-42 mm 
-300 to 750 mm 
> lxlO^''cm^^ 
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bulk currents have increased. The spread of the measured currents is partly due to the 
difference in the sensor positions relative to the interaction point, but is dominated by 
variations in the sensor temperatures. The occasional dips in the currents are related to 
short annealing periods in which the cooling systems were switched off. 

The measured currents contain contributions from two dominant sources, generically 
referred to as bulk and surface currents. The bulk currents vary exponentially with 
temperature and have a precisely predicted relationship with fluence. Surface currents 
arise due to irregularities introduced during sensor production such as process errors, 
scratches, guard rings and non-uniformities in the cut edges. Some of these contributions 
may have an exponential dependence on temperature [SJ, however the surface currents 
measured in VELO sensors are predominantly characterised by an Ohmic increase in 
current with bias voltage. In the majority of sensors the Ohmic surface current is seen to 
anneal with particle fluence. An analysis of VELO sensor currents as a function of bias 
voltage (IV scan) is described in detail in ref. 

For the VELO sensors the pre-irradiation exponential contribution is very small and is 
assumed to consist of a mixture of bulk and surface currents. The relative contribution of 
bulk and surface current is identified by measuring the current as a function of temperature 
(IT scan), as shown in figure [s] This method is described in detail in ref. [5]. 




2011 2011 2011 2011 2011 2011 

Figure 2: Currents measured in the VELO for each sensor as a function of time (bottom). The 
luminosity delivered to LHCb and the average sensor temperature is shown over the same time 
scale (middle and top). Increases in the dehvered luminosity are matched by increases in the 
sensor currents. The mean measured current agrees well with the prediction from simulation, 
which is described in Sect. 2. 4. The mean measured value excludes sensors that are surface- 
current dominated. 
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2.2 The effective band gap 



In the temperature range of interest, the bulk current is expected to scale according to, 

/(T)ocr^exp(-A,). (1) 

where T is the temperature in Kelvin and k is the Boltzmann constant (8.62 x 10^^ eV). 
The constant Eg is related to the bandgap energy and an effective value of 1.21 eV is 
assumed for the temperature range of interest [6] . Using this function the sensor currents 
can be normalised to a common temperature. In addition, the current variation as a 
function of temperature can be fitted to measure Eg in the VELO sensors. A summary 
of the results obtained from several IT scans is presented in table [2j The statistical 
uncertainty quoted is the width of a Gaussian fitted to the distribution of the measured 
values from all sensors. The largest source of systematic uncertainty is related to the 
accuracy with which sensor temperatures are known, and the choice of the temperature 
fitting range . 

The weighted average of the measured values is Eg = 1.16 ± 0.03 ± 0.04 eV, which is 
statistically compatible with the expected value of 1.21 eV from the literature [H]. Recent 
studies [0] have shown that discrepancies can occur due to sensor self-heating, in which case 
Eg is measured to be systematically high. In addition, dependencies of Eg on the sensor 
bias voltage have been observed. However, the VELO sensors used for these measurements 
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Figure 3: The current versus temperature for two VELO sensors operated at the nominal 150 V 
bias. Before irradiation there is a small exponential component for the sensor shown in figure a), 
while that shown in figure b) is dominated by the non-temperature dependent surface current 
contribution. After irradiation, a large exponential component is seen for both sensors. These 
sensors are now bulk current dominated and the surface current has largely been annealed. 
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Table 2: The effective band gap, Eg, measured following various amounts of delivered luminosity. 
The first uncertainty is statistical and the second is systematic. 



Delivered luminosity 


Bias voltage 




[fb-M 


[V] 


[eV] 


0.48 


100 


1.17 ±0.07 ±0.04 


0.48 


150 


1.18 ±0.05 ±0.04 


0.82 


150 


1.14 ±0.06 ±0.04 


1.20 


150 


1.15 ±0.04 ±0.04 



were cooled and sufficiently biased at 150 V, such that these effects should not significantly 
infiuence the result. This is supported by the consistency of the cross-check measurement 
made at 100 V. 

2.3 Fluence determination 

With a good understanding of the annealing conditions, changes in current can be related 
to the particle fiuence incident on a sensor. The expected change in leakage current at 
room temperature is given by the relation, 

AI = a<pVs^, (2) 

where a is the annealing parameter in units Acm~^, (j) is the fluence in units of number 
of particles per cm^, and Vsi is the silicon volume in cm^. The annealing parameter is 
a logarithmic function of time, and also depends on temperature. It has been shown to 
follow an Arrhenius relation, 

a OC exp 

where Ea is the activation energy for which a value of 1.33 ± 0.07 eV has been derived 
from a flt [7j. In order to estimate the appropriate value of a for our data, we proceed as 
follows. The sensor temperatures are recorded for each minute of operation and corrected 
to an equivalent time at 21 °C. The delivered luminosity, C, is folded with this equivalent 
time to produce an effective value of acp to be used in eq.[2] This procedure [3] yields 
a value for aC of 4.8 x 10"^^ Afb~^ cm~^, corresponding to an effective value for a of 
6 X 10-^^Acm-^ 

Fluences simulated with GEANT4 in the LHCb detector are folded with particle dam- 
age factors using the NIEL scaling hypothesis [H]. The damage caused by high energy 
particles in the bulk material is parameterised as proportional to a displacement damage 
cross section, D. A displacement damage cross section of 95 MeVmb for 1 MeV neutrons, 
Dn(lMeV), is assumed and each particle, i, is assigned an energy-dependent hardness 
factor, ki, 

UE) = D,{E)/D^{lMeV), (4) 
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which can be used to estimate the damage in the material. The total damage is expressed 
as a multiple of the expected damage of a 1 MeV neutron (referred to as 1 MeV neu- 
tron equivalent, or IMeVueq). This technique has been shown to be highly effective for 
describing the evolution of the leakage current. For this analysis, every particle in the 
simulation that is incident on a VELO sensor is assigned a displacement damage cross- 
section, or cross-section per silicon atom, using as reference a set of tabulated values from 
ref. 0. 

The fluence varies with the geometric location of the sensors in the VELO, as shown 
for simulated events in figure |4| It decreases with radial distance from the beam with 
an approximate dependence. The position along the beam-pipe (2;-direction) and 

path length of the particle in the silicon sensor were both found to significantly affect 
the fluence incident on a sensor. The predicted fluence is assigned an uncertainty of 
~8% which is dominated by uncertainties due to particles with no available displacement 
damage data. The average measured current increase for the sensors at the operational 
temperature of approximately — 7°C (and including several annealing periods) is ~18 [lA 
per fb~^. 

2.4 Predicted currents 

With the relationship between luminosity and damaging fluence established, the expected 
change in leakage current due to radiation damage can be predicted. The integrated lu- 
minosity is taken from the LHCb online measurement. Predictions are shown to agree 
with the measured currents, as shown in flgure [5] The predictions have an associated un- 
certainty of approximately 10%, estimated by adding the uncertainties for the integrated 
luminosity (5%), the annealing factor (3%) and the damaging fluence prediction (8%) 
in quadrature. The predicted leakage current is on average within 5% of the measured 
current. 

3 Depletion voltage studies 

The depletion voltage, deflned as the reverse bias voltage required to fully deplete a sen- 
sor, has been monitored as a function of the delivered luminosity. For fluences delivered 
to the VELO within the flrst few years of operation, the change in depletion voltage 
for an n^-on-n type sensor should be accurately described by the Hamburg model [TO] . 
The effective doping of the n-bulk changes over time due to radiation-induced defects. 
Dominant mechanisms are expected to be the inactivation of phosphorous dopants in 
combination with the introduction of acceptors. This leads to an initial decrease in the 
depletion voltage of the n^-on-n sensors to a value close to V. The n-type bulk inverts 
and becomes p-type, after which further irradiation leads to an increase in depletion volt- 
age. Eventually, the bias voltage required to obtain sufficient charge from a sensor will 
cause electrical breakdown in the silicon or exceed the 500 V hardware limit, thus limiting 
the useful lifetime of the silicon detector. For an oxygenated n^-on-p type sensor irradi- 
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Figure 4: a) The fluence from Ifb^^ of integrated luminosity versus radius for two VELO 
sensors, as seen in simulated proton-proton collisions at a 7TeV centre-of-mass energy. Top b) 
In each sensor, the fluence as a function of radius is fitted with the function Ar'^. The fitted 
exponent, k, is shown as a function of the sensor z-coordinate, where z is the distance along the 
beam-axis that a sensor is from the interaction region. The distribution of the fluence across the 
sensor is seen to become flatter with increasing distance from the interaction region. Bottom 
b) The fluence at the innermost radius of the sensor against the sensor z-coordinate. 

ated with charged hadrons there are expected to be competing mechanisms, with acceptor 
introduction partially compensated by initial oxygen-induced acceptor removal pTl IT2] . 

Following manufacture the depletion voltage of each VELO sensor was measured by 
comparing the capacitance (C) to the bias voltage (V) [13]. It is not possible to implement 
this technique after VELO installation and so alternative methods are used to extract 
information on the depletion voltage. This section presents results from two such studies: 
Charge Collection Efficiency and noise scan studies. 

3.1 Charge Collection Efficiency 

The amount of charge collected by an under-depleted silicon strip increases with the bias 
voltage applied. When the sensor is fully depleted, any further increase in bias voltage 
will not increase the amount of charge collected (given a sufficient signal collection time). 
The relationship between the Charge Collection Efficiency (CCE) and the applied bias 
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Figure 5: The leakage current against sensor z-coordinate after 1.20 fb~ of integrated luminosity, 
normalised to °C. The data is in agreement with predictions, represented by the shaded region. 
The two VELO halves are referred to as the A and C sides of the VELO. 



voltage has been exploited to measure a property of the sensor analogous to the depletion 
voltage, referred to as the Effective Depletion Voltage (EDV). 

3.1.1 Effective Depletion Voltage determination 

The nominal operational voltage of the VELO sensors is 150 V. For the CCE analysis, 
collision data is recorded with every fifth module operated at a voltage ranging between 
and 150 V. The remaining modules are maintained at 150 V. Sensors with variable voltage 
are referred to as test sensors. The test sensors are removed from the reconstruction 
algorithms such that only hits from the 150 V operated sensors are used to reconstruct 
particle tracks. A track is extrapolated to a coordinate on the test sensor and the set 
of five strips nearest to this coordinate are searched for deposited charge. This provides 
unbiased information on the amount of charge deposited by the particle as a function 
of bias voltage. At each bias voltage the pedestal-subtracted ADC distribution is fitted 
using a gaussian convoluted with a Landau function. This is used to determine the Most 
Probable Value (MPV) of the ADC distribution. At large bias voltages the MPV of the 
ADC distribution reaches a plateau. The EDV is defined as the voltage at which the MPV 
of a sensor is equal to 80% of the plateau ADC value, as shown in figure [6j The threshold 
of 80% was chosen as it gives closest agreement with depletion voltages determined from 
pre- irradiation CV measurements [13] . For unirradiated sensors, the difference in the 
values obtained using the CV method and the EDV method is less than 10 V. 

Differences between the depletion voltage and EDV are expected near to sensor type- 
inversion. The depletion voltage is expected to decrease to a value of approximately V, 
whereas the minimum EDV is dictated by the smallest potential difference required to 
collect charge from the silicon strips, which in turn depends on the shaping time of the 
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Figure 6: a) The pedestal subtracted ADC distributions for an R-type sensor at three example 
bias voltages, b) The MPV of the fit to the ADC distribution vs. bias voltage. The dashed 
lines represent the ADC that is 80% of the plateau value, and the corresponding EDV. 



electronics. When operated with a 25 ns signal shaping time, the smallest EDV that can 
be measured in the VELO sensors is approximately 20 V. The EDV is therefore not an 
accurate estimate of the depletion voltage below this value. 



3.1.2 Bulk radiation damage 

Between April 2010 and October 2011 five dedicated CCE scans were taken, corresponding 
to delivered luminosities of 0, 0.04, 0.43, 0.80 and 1.22 fb~^. As the fiuence delivered to 
the sensors varies significantly with sensor radius, each sensor is divided into 5 radial 
regions such that the fiuence does not change by more than a factor of two across a 
region. The change in EDV with irradiation for a particular n^-on-n type sensor is 
shown in figure |7(a)[ Initially the EDV is found to decrease with fiuence across all radial 
regions, as predicted by the Hamburg model. The rate of decrease is greater in the inner 
radial regions of the sensor, consistent with expectations that these regions are exposed 
to higher fiuence. The innermost region undergoes an increase in EDV between 0.80 and 
1.22 fb~^ of delivered luminosity, indicating that this part of the sensor has undergone 
type inversion. The n~^-on-p type sensors exhibit a decrease in EDV with initial fiuence. 



as shown in figure 7(b) This is understood to be caused by oxygen induced removal of 
boron interstitial acceptor sites, an effect that has been previously observed [HI [12] . 
The global change in EDV is determined by combining the data from many of the 



VELO sensors with the predicted fiuence (see section 2.3), as shown in figure [8j Sensors 
are divided into categories based on their initial EDVs. The irradiation-induced change 
in the depletion voltage of n"'"-on-n type sensors is modeled as a function of time, tem- 
perature and fiuence by the Hamburg Model. It has three components: a short term 
annealing component, a stable damage component and a reverse annealing component. 



9 




30 35 40 45 50 
Sensor radius [ mm ] 



30 35 40 45 50 
Sensor radius [ mm ] 



(a) 



(b) 



Figure 7: a) The EDV against sensor radius for an n^-on-n type sensor for each of the CCE 
scans. The dashed hne shows the mean EDV across all radius regions prior to sensor irradiation, 
where some Ofb~^ data points are not present due to low statistics, b) A similar plot for the 
n~^-on-p, (p-type sensor. The minimum EDV is ~40V, which is significantly higher than the 
minimum at ~20V observed for the n"'"-on-n type sensor. 



Taking into account the LHCb luminosity measurements and VELO sensor temperature 
readings, the Hamburg model predictions can be compared to data, as shown by the 
overlaid curves in figure [8j Good agreement is found for low fiuences, and for higher 
fluences after type inversion. It is assumed that the sensors type invert at a fiuence 
near to the EDV minimum. For all n"*'-on-n type sensors this occurs at approximately 
the same fiuence of (10 — 15) x lO^^lMeVuoq. The behaviour after inversion is found 
to be independent of the initial EDV of the sensor, with an approximately linear in- 
crease in EDV with further fiuence. A linear fit to the data gives a voltage increase 
with fiuence of (1.35 ± 0.25) x 10~^^ V/1 MeVueq. For the n~^-on-p type, the initial de- 
crease in EDV occurred up to a fiuence of approximately 2 x 10^^ IMeVugq. After this 
the EDV has increased with further fiuence. The rate of increase is measured to be 
(1.43 it 0.16) X 10~^^ V/1 MeVueq, similar to that of the type inverted n^-on-n type sen- 
sors. 

The EDV of the n^-on-p type sensors begins to increase having received significantly 
less fiuence than the n"'"-on-n type sensors. If the comparable rate of EDV increase is 
maintained with further fiuence then the n^-on-p type sensors will reach an EDV of 500 V, 
the hardware limit of the VELO system, after receiving approximately 35 x 10^^ 1 MeV Ugq 
less fiuence than the n^-on-n type sensors. It is expected that the n"*"-on-n sensors will 
reach the 500 V limit following a fiuence of approximately 380 x 10^^ 1 MeVueq. 

The amount of charge collected is expected to change with fiuence due to radiation 
induced changes to the silicon. For 0-type sensors the MPV has decreased by approxi- 
mately 4% in the most irradiated regions, having received a fiuence of 40 x 10^^ 1 MeV Ugq. 



10 




Figure 8: The EDV against fluence for VELO sensors of various initial EDV. The EDV from 
data is compared to depletion voltages predicted by the Hamburg model, with good agreement 
observed prior to, and after sensor type-inversion. 



An even larger reduction of approximately 8% is found in the inner regions of the R-type 
sensors, having received a comparable fluence. This is due to a charge loss mechanism 
related to the second metal layer of the R-type sensors, which is described in detail in 
section |4j The outer regions of the sensor are most significantly affected, with decreases 
of approximately 12% observed following a fluence of just 2 x 10^^ IMeVueq. 

3.2 Noise scans 



The CCE scan data described in section 3.1 requires proton beams, and so is collected at 



the expense of physics data. A second method has been developed to monitor radiation 
damage, using the relationship between the intrinsic electronic noise of the pre-amplifier 
and the capacitance of the sensor. Data scans for this study can be collected regularly as 
proton collisions are not required. 

In undepleted silicon, several sources of input capacitance are identified, the most 
dominant of which is the inter-strip impedance. For n^-on-n sensors before type inversion, 
the depletion region grows with increasing voltage from the backplane (the opposite side 
to the strips). When the sensor is fully depleted the space-charge reaches the strips 
and the inter-strip resistance increases by several orders of magnitude, resulting in a 
decrease in sensor noise [HI [15]. For n'*'-on-n type sensors following type inversion and 
n'^-on-p type sensors, the depletion region grows from the strip side of the silicon. In this 
situation the strips are immediately isolated at the application of a bias voltage and the 
relationship between noise and voltage cannot be exploited to extract information related 
to the depletion voltage. 
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The intrinsic noise in VELO sensors is determined by subtracting the mean ADC 
value (or pedestal) and a common-mode noise term. Figurejo] shows the inverse of the 
intrinsic noise as a function of voltage for an n'^-on-n and n~^-on-p type sensor. For the 
n'^-on-p before and after irradiation, and n"*"-on-n after irradiation (having type inverted) 
the distribution is fiat, thus little information related to the sensor depletion voltages can 
be extracted. For the n^-on-n type sensor prior to type inversion, an increase in voltage 
results in a decrease in noise until a plateau is reached when the sensor is fully depleted. 

The noise scan data can be used to identify whether an n'*'-on-n type sensor has 
undergone type inversion. Only R-type sensors are investigated as the strip orientation 
allows the identification of strips that have been subject to a specific fiuence. Following a 
delivered luminosity of approximately 0.80 fb~^, 40 of the sensors are identified as having 
type-inverted in the first radial region (8—11 mm), while in the second region (11 — 16 mm) 
21 sensors are identified. Similar information can be extracted from the CCE data, with 
a sensor region defined as type-inverted when the measured EDV has reached a minimum 
and subsequently begun to increase. Following the same luminosity, the CCE method 
identified 21 and 5 type-inverted sensors in the first and second radial regions, considerably 
fewer than the noise method. This discrepancy is understood by examination of figure |8} 
in which the minimum of the Hamburg model prediction and the point at which the EDVs 
begin to increase are separated by a fiuence of approximately 10 x 10^^ IMeVueq. The 
noise scan method is not subject to the same fiuence lag. Following a delivered luminosity 
of 1.2 fb~^, the CCE method identifies 39 and 21 sensors in the two radial regions. This 
is in good agreement with the noise method, with the same 39 and 21 sensors identified 
by each method. 
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Figure 9: The inverse of the sensor noise against bias voltage for a particular n"'"-on-n and 
n^-on-p type sensor, for two values of integrated luminosities. The C-V scan measured initial 
depletion voltages are shown by the dashed lines. 
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4 Charge loss to the second metal layer 



All physics analyses at LHCb rely on efficient track reconstruction using clusters from the 
VELO sensors. A cluster is defined as one or several adjacent strips with charge above a 
particular threshold. The data samples described in section 3.1 are also used to measure 
the Cluster Finding Efficiency (CFE), by looking for the presence of a cluster at the track 
intercept on the test sensor. Before irradiation the mean CFE of the VELO sensors was 
greater than 99% [16]. After irradiation the CFE in many sensors decreased significantly, 
as shown in figure 10 The inefficiency is particularly prevalent at large sensor radii and 
for high bias voltages. The decrease in CFE does not appear to be proportional to the 
delivered luminosity, but instead exhibits a rapid drop between the 0.04 and 0.43 fb~^ 
data scans. 

To determine the source of this CFE decrease, a large sample of regular LHCb physics 
data has been used to measure the CFE for small spatial regions on a sensor at the nominal 
150 V bias. The result of this is shown in figure [TH displayed beside a diagram illustrating 
the layout of the second metal layer readout lines. There is a clear correspondence between 
the two figures, with high CFE measured in regions that are devoid of second metal layer 
lines. The relative orientation of the strips and routing lines in R-type sensors is shown in 



figure 12 In addition, a schematic cross-section of an R-type sensor is shown in figure 13 



The routing lines from inner strips are seen to pass approximately perpendicularly over 
the outer strip implants. 

Using the precision tracking of the VELO, it is possible to investigate the CFE loss as 
a function of the distance between a track intercept with a sensor, and the nearest strip 



and routing line. This is shown in figure 14 The CFE is improved for track intercepts 



that are near to the strip implants. Conversely, the CFE is reduced when a track intercept 
is far from a strip and near to a routing line. Similar effects have been observed in other 
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Figure 10: The CFE for an R-type n~^-on-n sensor as a function of sensor radius for a) different 
amounts of delivered luminosity and b) several different bias voltages. 
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Figure 11: a) The layout of the second metal layer routing lines on an R-type sensor. The 
darker regions represent the presence of routing lines, and the lighter regions their absence, b) 
The CFE shown in small spatial regions of an R-type sensor. 



experiments (TTJ [18] . 

The source of the CFE loss is hypothesised in terms of charge induction on the second 
metal layer. Prior to irradiation, ionised electrons will drift along the field lines, most of 
which terminate at the implants. Hence the majority of the signal will be induced on 
the implants, which are strongly capacitively coupled to the readout strips. The drifting 
charge is expected to be collected well within the ~20 ns readout period of the electronics, 
and no signal is expected on neighbouring electrodes (with the exception of capacitive 
coupling and cross-talk effects, which are measured to be low). However, irradiation may 
cause modifications to the field line structure, such that not all field lines terminate on 
the implants. In addition, there may be charge trapping effects which delay the drift of 
charge, resulting in charge sampling before the electrons have reached the implants. In 




Figure 12: A photograph of the innermost region of an R-type sensor. Strips run from the 
bottom-left to the top-right. Each strip is connected to a routing line orientated perpendicularly 
to the strip. 
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Figure 13: A schematic cross-section of a portion of an R-type sensor, showing the relative 
position of the two metal layers used to carry the readout signals in n^-on-n type sensors. The 
implants and strips (into the page) run perpendicularly to the routing lines (left to right). 
For clarity the routing line of just one strip is shown. 
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Figure 14: The CFE as a function of the distance between the particle intercept and the nearest 
routing line, for several bins of the distance between the particle intercept and closest strip edge. 



both of these situations there will be a net induced charge on nearby electrodes, such as 
the second metal layer routing lines. In figure 10(b)| the CFE was seen to worsen with 
increasing bias voltage. This appears to disfavour the contribution due to trapping, as an 
increase in bias voltage should result in faster collection times. However, the bias voltage 
may also affect the field-line structure. In reality, it is likely that the charge loss to the 
second metal layer is due to several competing mechanisms. 



The CFE loss also exhibits a significant radial dependence, as was shown by figure 10 



This can be understood by considering two competing mechanisms. The implant strip 
width and the fractional area covered by the strips increases with radius, resulting in 
reduced charge loss, due to greater strip shielding. However, the fractional area covered 
by the second metal layer also increases with radius, due to the greater density of lines, 
increasing the amount of pickup. The latter effect is dominant, hence the overall charge 
loss is greater at large sensor radii. 

In addition to lowering the clustering efficiency, charge induced on a routing line can 
introduce a noise cluster. The cluster ADC distribution from R-type sensors has a peak 
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Figure 15: For R-type sensors: a) The ADC spectrum of all clusters seen for three different 
integrated luminosities. The limit at 10 ADC counts is imposed by the clustering thresholds, 
b) The fraction of reconstructed clusters that are induced on routing lines as a function of the 
distance to the nearest routing line and strip. This is determined from the number of track 
intercepts for which the inner strip associated to the nearest routing line has a 1 strip cluster 
with less than 35 ADC counts. 



associated to these low ADC noise clusters that has grown with fluence, as shown in 



figure 15(a), These noise clusters are predominantly single strip clusters located at small 
radius regions of R-type sensors. The fraction of the induced noise clusters increases 
when tracks traverse a sensor near to a routing line and far from a strip, as shown in 



figure 15(b) 



The CFE decrease is not observed in 0-type sensors as the routing lines from inner 
strips were intentionally placed directly above the outer strips to minimise pick-up. This 
is made possible by the 0-type sensors strip orientation. The CFE loss could be partially 
recovered by lowering the cluster reconstruction thresholds. However this comes at the 
expense of a worse signal to background ratio which leads to higher rates of fake tracks 
reconstructed from noise induced clusters. 

The R-type n~^-on-p and n^-on-n sensors have a similar CFE dependence on the strip 



and routing line distance, as shown in figure 16(a) Figure 16(b) shows the MPV of the 



collected charge distribution as a function of bias voltage for an n+-on-n and an n~^-on-p 
type sensor. At 150 V, the MPV of the n"'"-on-n and n^-on-p are approximately equal, 
both before and after irradiation. For the n"'"-on-n type sensor post irradiation, the MPV 
reaches a maximum at around 60 V after which it is observed to decrease with increasing 
bias voltage. This decrease in MPV leads to a reduced CFE and is associated with the 
second metal layer effect. Therefore less charge is lost to the second metal layer when 
operating the n^-on-n sensor at a lower than nominal voltage. 

The n~^-on-p type sensor does not exhibit this same charge collection loss (and resulting 
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Figure 16: a) The CFE as a function of the distance between the particle intercept and the 
nearest routing line and strip edge, compared for an n"''-on-n and n~^-on-p sensor, b) The MPV 
vs. bias voltage for an n+-on-n and n~^-on-p sensor. 
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Figure 17: The MPV against bias voltage for an R-type sensor at three values of delivered 
luminosity. The sensor region has been identified as having type inverted in the 1.20 fb~^ scan, 
in which the MPV dependence on voltage is no longer present. 



CFE decrease) dependence on voltage. This may be due to the depletion region in the 
n^-on-jo type sensor growing from the strip side of the silicon instead of from the sensor 
backplane. This is supported by the observation that after type inversion, charge loss in 



n"*"-on-n type sensors no longer depends on the bias voltage, as shown in figure 17 
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5 Summary 



The effects of radiation damage have been observed in all of the LHCb VELO sensors. 
At the operational sensor temperature of approximately — 7°C, the average rate of sensor 
current increase is measured to be 18 |J.A per fb~^, in agreement with expectations. The 
silicon effective bandgap has been determined using current versus temperature scan data 
collected at various levels of radiation, with an average value of Eg — 1.16±0.03±0.04eV 
found. 

Analysis of the Charge Collection Efficiency has proven an effective method for track- 
ing the evolution of the sensor depletion voltages with fluence. Measurements of the 
Effective Depletion Voltage after sensor type inversion have been shown to agree well 
with the Hamburg model predictions. A method relating the sensor noise dependence 
on bias voltage was used to identify type-inverted sensors, and was found to be in good 
agreement with the Charge Collection Efficiency method. 

A significant decrease in the Cluster Finding Efficiency in R-type sensors due to the 
second metal layer has been observed following a relatively small amount of particle 
fluence. In the worst affected sensor regions the Cluster Finding Efficiency decreased by 
over 5%. Despite this relatively large localised inefficiency, studies of the VELO tracking 
efficiencies show no degradation associated to this effect, within the errors of ±0.3%. For 
the n^-on-n type sensors before type inversion the magnitude of the charge loss is found 
to increase with sensor bias voltage. For type-inverted n+-on-n type sensors and n^-on-p 
type sensors, a voltage dependence is not observed. Radiation induced modifications to 
the field line structure and charge trapping effects are thought to be possible sources of 
this charge loss effect. 

The two n~^-on-p sensors have been studied in detail, providing valuable information 
for upgraded detector designs. If the Effective Depletion Voltage were to continue increas- 
ing at the currently observed rate with further irradiation, then the n'^-on-p would reach 
the 500 V hardware limit having received approximately 35 x 10^^ IMeVueq less fluence 
than an equivalent n"'"-on-n type sensor. This corresponds to approximately 1 fb~ of 
integrated luminosity in the highest particle fluence region of the VELO. 
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